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ABSTRACT 


Displacement tests have been conducted on an un- 
consolidated Ottawa silica sand pack to study the effect 
of viscous fingering on the recovery of oil at water 
breakthrough. 

Tests were performed with three oils using dis- 
tilled water as the displacing fluid. The oil-water vis- 
cosity ratios were 15.5, 34.2 and 111.4 respectively. 

The recovery data were correlated using the dimensionless 
group I, which represents the ratio of viscous to cap- 
illary forces. 

The results indicate that the value of I at which 
the effect of viscous fingers manifests itself as a 
reduction of the breakthrough recovery may be estimated 


from: 
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if the constant, C, is taken as 30 for each system. 
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1. INTRODUCTION 


Waterflooding has by far been the most useful second- 
ary recovery technique developed by the oil industry to 
increase the recoveries from depleting oil reservoirs 
around the world. Extensive research and past experience 
have provided methods to predict production performance 
within a reasonable confidence level. There is however, 

a constant need for improvements as the industry takes on 
Marginal projects with a rapid increase in capital costs. 

It is therefore important to study some of the uncer- 
tainties which still exist regarding the effect of certain 
variables. Rate is probably the most important of these, 
Since it is fairly easy to control and has a huge impact 
on the economics of a project. 

The interest in recovery of heavy high viscosity 
crudes is rapidly gaining momentum. The literature deal- 
ing with immiscible displacements at such unfavourable 
viscosity ratios contains obvious inconsistencies. An 
especially controversial subject has been the effect of 
viscous fingering on the displacement efficiency. Al- 
though it may not be a significant phenomenon in an actual 
reservoir, fingering must be taken into account when 
evaluating data from otherwise scaled laboratory models. 


This study was conducted to obtain a better 
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understanding of the viscous fingering phenomenon and its 
effect on the recovery of oil from an unconsolidated 
porous medium. Three different oils with a wide range in 
viscosity were employed; the rate of displacement was 
varied from 2.5 to 3520 cc/hr and recoveries at break- 
through and Suh eSquent production behaviour were recorded 


and correlated using a classical scaling parameter. 
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2. LITERATURE REVIEW 
2.1 Viscous Fingering 


Fingering of the advancing water into the less mobile 
Oil was early recognized as a problem in the oil industry. 
However, it was mainly attributed to gravity effects, 
permeability stratifications and other gross inhomogen- 
eities. 

Miller? reported in 1941: "Visual observations of 
the oil-filled test sand column during the water encroach- 
ment experiments indicate conclusively that different 
volumes of oil are recovered from different unit volumes 
of sand as a result of 'fingering' of the water and bypas- 
Sing of 011". However, he made no special point of this 
in his discussion of the data. 

Engelberts and Hi inkenbera.. in 1951, were probably 
the first to point out the important role the instability 
of the transition zone played when a fluid is displaced by 
one which has a lower viscosity. They called the phenom- 
enon "viscous fingering". Their results showed a consis- 
tent decrease in breakthrough recoveries with increasing 
viscosity ratios. They also observed a decrease in the 
breakthrough recovery with increasing rates (Uo/ Hy = 24) 
for both oil-wet and water-wet media. 


Van Meurs and van der Poel? provided a simple 
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mathematical description of a water drive involving vis- 


cous fingering. Their equations can be used to predict 
the production performance in a field with regularly 
spaced wells. The data required are; the viscosity ratio, 
the residual oil saturation and the minimum water satura- 
tion. A water-cut value must be selected to complete the 
evaluation of a specific project. They obtained satis- 
factory agreement with production data from an actual 
field. The theory is presented in detail in Chapter 3. 
Chuoke et Ss presented a theoretical description 
of the instability of fluid displacements in porous media. 
They derived an equation representing the wavelength of 
maximum instability, or the most probable peak-to-peak 
distance, for a given system when fingering occurs. 
Assuming neutral wettability and horizontal displacement 


(i.e. gravity effects can be ignored) the equation becomes: 


Meee Tame mY 
where 

han = the most probable finger spacing (cm) 

C =a dimensionless constant 

yY = interfacial tension (dynes/cm) 

K = absolute permeability (cm?) 

V = total flow rate per unit cross-sectional 


area or volumetric velocity (cm/sec) 
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Wea viscosity of oil (poise) 


v1 viscosity of water (poise) 


Ww 
Instability will occur if the volumetric velocity exceeds 
a critical value for a specific rock-fluid system, pro- 
vided the perturbation contains wavelengths greater than 
a critical wavelength, hoe The relationship between the 


two wavelengths is: 
(2) 


where 
eo = the critical wavelength or finger distance (cm). 
An expression for We as a fraction of the tube diameter is 


obtained by a combination of equations 1 and 2: 


r 
ares 3.V ee )D* | ‘ 
Oo "Ww 

where 

D = the diameter of the tube (cm). 
From the above definition of hoe it follows that the oil 
water interface is stable for any rate if the critical 
finger distance exceeds the diameter of the tube. Equation 
3 would have been a powerful tool if the constant C could 
be easily predicted. Chuoke et al. obtained a reasonable 


fit to their observed finger spacing if C in equation 1 


was taken as 30. That was for a neutral wettability 
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system containing no initial water. Although their exper- 

imental data were insufficient to calculate a value, they 

indicated that C was larger for water-wet media containing 

connate water. The value of C must be established experi- 

mentally for each system. Let us assume that viscous 

fingers will be formed as soon as he becomes equal to D. 

If this shows up in the recovery data as a noticeable 

decrease, then this point can be used to determine the con- 

stant in equation 3 from: 

3.V (uo=n,,)D* 
Oo 'W 

YK 


(4) 


Chuoke et al. stressed the fact that larger values of their 
constant lead to increased critical wavelengths. The lateral 
extent of a laboratory model or core sample may then easily 
be exceeded. They suggested that this might explain why 
viscous fingering effects may not reveal itself even though 
the system satisfies all other criteria for instability. 
Their visual observations showed that increased viscosity 
ratios and decreased interfacial tensions both favoured 

the formation of smaller and more numerous fingers. 

The main limitation on the applicability of their 
theory is the assumption that the Muskat model of dis- 
placement, where oil and water flow in separate macro- 
scopic regions, is valid for an actual reservoir contain- 


ing an initial water saturation. 
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de Haan” provided experimental support for Chuoke 
et al.'s* theory. He was able to predict the rate at 
which fingers started to form using their equations. The 
recoveries at low rates were high, decreasing in an inter- 
mediate range and levelling out at high rates. He con- 
cluded that at low rates, water tended to fill up the 
whole cross-section, and an efficient frontal displacement 
resulted. In the intermediate range, fingers were formed, 
and as a result the recoveries dropped. At high rates, 
numerous fingers were observed, and the recovery became 
insensitive to further rate increases. 

de Haan? explained that viscous fingering effects 
are unlikely to be pronounced in oil-wet systems due to 
the distribution of the fluids. The constrictions are 
filled with oil in an idealized oil-wet system, and the 
flood fronts in adjacent channels will anenerore move 
quite independently. Fingers of macroscopic dimensions 
will only be formed when a more or less coherent oil- 


water interface exists, and this condition is clearly not 


met in an oil-wet system. 


Cir dealt with the instability problem 


Scheidegger 
in two theoretical papers. According to his analyses 
fingering should be independent of the displacement vel- 


Ocity. He cited the work of Blackwell et Bie to support 


this controversial conclusion. These authors obtained 
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breakthrough recoveries that were independent of rate with 
a Mobility ratio of°93."*An explanation” of this ‘surprising 
result might be the fact that the displacement rates used, 
1.0 - 100 ft/day, placed all their results in the insensi- 
tive range described by de Haan>, 

Outmans” included some of the nonlinear terms prev- 
iously neglected in the analyses of the instability 
equations. He showed that viscous fingering then could 
become independent of rate only if both gravity and inter- 
facial tension effects were negligible. 

Rachford?° was the first to use the Buckley-Leverett 
flow model, in which oil and water flow simultaneously, in 
the analyses of the stability of the interface. He includ- 
er the effects of connate water and the transition zone, 
both normally present when water-wet systems are flooded. 
The transition zone tended to insulate incipient fingers 
from the high-mobility water and thus reduced the finger- 
ing effect. This is not accounted for in.the parallel 


plate theory. 


My 10 


Perkins and Johnston confirmed Rachford's 
findings experimentally. The fingering behaviour differed 
in the presence of connate water. Numerous small fingers 
developed near the entrance, but they soon broke up and 
formed a graded saturation zone. No sharp front could be 


followed. Perkins and Johnston suggested that the dampening 
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of the fingers was caused by crossflow of the phases, 
transverse to the direction of the gross fluid movement. 
Their experiments supported this explanation, and they 
concluded that immiscible dispersion ought to be incorp- 
Orated in the analyses of the stability to adequately 
describe the viscous fingering phenomenon in natural 
systems. 

Hagoort?? highlighted the effect Hee the choice of 
flow-model (i.e. Muskat or Buckley-Leverett) has on the 
stability criteria. He showed that the Buckley-Leverett 
approach, previously only attempted by Rachfora’®, yields 
different criteria for the initiation of fingers than those 
derived from a Muskat-model. He found that the interface 
became unstable if the so-called shock mobility ratio 
was greater than 1, provided the wavelength of the instab- 
ilities was smaller than the canal width. . The shock mob- 


ility ratio is defined as: 


nen the mobility of fluids behind the shock front 
S the mobility of fluids ahead of the shock front 


Or 
M = Koow Kro (So) 7, : Kor Kew (Sy 
‘ Kecw "o 
where 
Ke = permeability to oil at the connate water 
saturation 
Kg (Sa) = normalized relative permeability to oil at 


the shock saturation 
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Sa = permeability to water at the residual oil 
saturation 
RUS) = normalized relative permeability to water 


at the shock saturation 


The shock mobility ratio is considerably lower than the 
end-point mobility ratio which determines the stability of 
a Muskat-type displacement. The higher the mobility ratio, 
the faster instabilities arise. The Buckley-Leverett 
approach is generally believed to describe immiscible 
displacements in water-wet media better than the Muskat- 
model. The stability of these displacements is therefore 
probably better than earlier results based on the Muskat 


flow model indicated. 


2.2- Scaling 


Numerous papers on the scaling of waterflood models 
have been published since Leverett et she ee introduced the 
concept in 1942. In the present work no effort was made 
to scale any particular reservoir. The problem was there- 
fore limited to the end-effects and the fingering phenom- 
enon. 

The dimensionless scaling group, I, defined as: 

LV Uy 
yV/K_ 
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11 
where 
L = length of the system _ (cm) 
was used by Engelberts and Klinkenberg” to correlate hor- 
izontal displacement experiments. It is simply a measure 
of the ratio of the viscous to the capillary forces. 


14 and 


de Haan”, uSing the data of Rapoport and Leas 
Kyte and Rapoport!, concluded that the recovery became 
independent of I when I was larger than 0.1. Substitution 


of normal field-values indicated that I was always in this 


insensitive range. 


The same scaling group has been used by Jones-Parra 


et al +o Collins’, scott !® and Kloepfer?? 


which facili- 
tated the comparison of results. 
Geertsma et eioees arrived at a more universal scaling 


number given by: 


L = Weos 9 VK¢ 


t Dee 
where 
8 = wetting angle of fluid interface (dimensionless) 
d = porosity (dimensionless) 


This is the reciprocal: of I, except that the effects of 
porosity and contact angle are included. The problems in- 
volved in establishing the contact angle makes the applic- 


ation of this scaling group difficult. 
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End-effects, as described by Kyte and Rapoport!>, 
can Obscure data obtained from unscaled laboratory models. 
Spontaneous localized imbibition or the inlet effect is 
minimized if rate and viscosity ratio are kept low. 

This effect could probably be eliminated if the inlet 
system was designed to assure uniform liquid contact at 
the sandface. The so-called outlet effect, caused by the 
discontinuity of the capillary pressure, becomes less 
severe when rates and viscosity ratios are increased. The 
relative importance of both effects decreases if longer 
systems are employed. 

Chuoke et ae indicated that scaling of viscous 
fingering should be straightforward in most cases. Prob- 
lems would arise however, if the wavelength of maximum 
instability was of the same magnitude as the model width. 

Outmans” indicated that it might prove impossible 
to accurately scale the fingering phenomenon, because to 
achieve similarity the small heterogeneities that initiate 
the fingers have to be scaled. 

Rachford~° concluded that no additional scaling 
requirements were necessary to deal with the instability 
problem. His model, based on the Buckley-Leverett dis- 
placement mechanism assuming a porous medium with a con- 
nate water saturation, should be more representative than 


the previous idealized parallel plate models. 
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Perkins and Johnston ~ showed that any experiment 
aimed at scaling immiscible fingering in a water-wet 
reservoir should be performed with an initial connate 
water saturation. However, no indication was given as to 


the need for matching the initial saturation of the proto- 


type. 
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3. THEORY 


van Meurs and van der Poel? presented a flow-model 
that accounted for viscous fingering in immiscible dis- 
placements. Their simplifying assumptions were based on 
pictures of scaled waterflood experiments in a transparent 
system. According to their observations the flow mechan- 


ism could be idealized in the following way. 


1. The displacement process is of a one dimensional 
character in the region where both oil and water 
are present (i.e. movement is parallel to the 
direction of flow). 

2. Oil flow outside the water fingers and water flow 
in the center of the fingers proceeds unhindered. 

3. In the edge zone of the water fingers both water 
in protrusions and oil in enclosed pockets are 
completely prevented from moving in the direction 
of flow. 

4. The fractions of both immobile oil and water in 
a cross section of the formation are independent 


of time and place. 


They considered a unit cross section of the formation 
large enough to contain sufficient fingers to be represent- 
ative of the total cross section, see Figure 1. This 


section was divided into three areas: 
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a) Area 1, the sum of all areas outside the water 
fingers where oil flows unhindered. 

b) Area 2, the sum of the center areas of the 
water fingers where water flows unhindered. 

c) Area 3, the sum of all areas containing edges 
and protrusions. In this area both the oil and 


the water are immobile. 


[] waTeR O1L 


FIGURE 1: IDEALIZED PICTURE OF WATER FINGERS 


The sizes of these areas can be characterized by 
means of oil and water saturations. 

All immobile water of Area 3 expressed as a fraction 
will be denoted Som" This means that eee is the average 


immobile water saturation for the unit cross section we are 
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looking at. In the same way ae is defined by averaging 
the immobile oil saturation. If the average water and oil 
saturations of the unit cross section are denoted Sy and 


So respectively, then the areas can be written: 


Area = "Se = 2S 
o) or 
Area 2=S -S§S 
Ww wm 
or since 
Sy rte Se isso 1 
WwW fe) 
Area l1=1- 85S - § 
or WwW 


LE Qs is defined as the flow rate of water per unit 
cross sectional area (dimension, L ee then the volumet- 


ric velocity of the water in Area 2 is defined by: 


W 5) =S 
w wm 


The flow in this area was assumed unhindered and Darcy's 


Law for a single phase can be used: 


dP 
Vena te K a (6) 
ay 
where 
sae = pressure in the water phase 
x = distance travelled in x - direction 


The same analyses for the oil phase yields: 
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and 
oP 
re) u ox 
fe) 
where 
Pe = pressure in the oil phase 


Although the authors did not mention it specifically, it 
should be noted that the use of the absolute permeability 
in equations 1 and 2 is justified by the assumption of 
unhindered flow. The oil present in Area 2 and the water 
present in Area 1 do not affect the ability of the porous 
medium to conduct fluids. Hence, the medium behaves as if 
it was completely saturated with a single fluid. Since 
macroscopic parts of the oil-saturated formation are 
bypassed by the penetrating water fingers, it is under- 
standable that the difference in pressure in the two 
phases, Pur Stays approximately constant as long as the 
finger edges are not interfering. This assumption was 


found to be true up to high water saturations. Consequent- 


ye 
P =P = P_ = constant 
fo: re) WwW 
and 
OP dP, 
re Boe (8) 
Ox ox 


from Equations 6, 8 and 9 it follows that: 


V.. = MV (10) 
WwW Oo 
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The total amount of fluid passing through the cross section 
is: 

CS Bee ® ah a 0 (Gr) 
and will be a function of time only since the fluids are 
assumed incompressible. 

An equation for the fractional flow of water is 
obtained by combining 5, 7, 10 and 12: 


M (S -S_) 
Pe i yh ORS tw WING. eal (13) 


eae es ae Gay(sh=siy sas 
where 


Boa b= So = 38 (14) 


dq. as 

= + ¢ ae = 0 (a5) 
or 

of 0S 

ot wan, (16) 
where 

Ws = a = cumulative water injected as a fraction 


of the pore volume 


L = Length of the formation 
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X = = the dimensionless distance along the flow 


ba 
L 
direction taken from the point of injection. 
ae is a function of sy only according to Equation 13. 
Equation 16 can therefore be written as: 
df Joe 0S. 


ee ee SO (ey) 
Ww oy 


Differentiation of Equation 9 gives: 


W MB 
Se ee eee (18) 
ds 2 

Ww [ (M-1) Eee ety, +B] 


The general solution of this differential equation is: 


ee w(S_) Bi W. (19) 


Li the initial water £ront (atst-= 0 or W, = 0) is 
assumed to be undisturbed and perpendicular to the flow 


direction then: 
w(S_) 9.0) (20) 


Substitution of Equations 18 and 20 in Equation 19 yields: 
MBW. 


X = ——___+_____—., (21) 
Wiieal Geist iy aceite 


and we can solve for the water saturation distribution 


in the region where both oil and water flow: 


Ww wm M-1 M-1 ax. (22) 


i 
{ { 
: val & enh 
WOL2 wrt wise) x iS iatacl ch cs aia 
VW ¢ ey fights) t : Taso 


- ‘jwLinnn= oF Sete 5 ts ee 


6 ngbfiaey 
a ~® 
vig 
in EV 
~ 
Mieranees sk Eee. SES AY: ep - 
ot: ; : 
y i . Y ‘ 2 
a Retvts. ed soie fabta Sh Ab 
h tpl fs ' “ah - . 7 
' here ae ha + (abl. 
: a) ond a% 
a - ; : 7 ’ = 


fh Oom iT mg okhy jpn aa 


7 Sine? 
of? o8¢ ad “rates or oer 144 lyn Bette se cai 
“ i | 7 : f 
a - oe 
\ i ’ a | 
\ . : 5 f ’ ‘ a 
iphloky @L or uf feipa. ri PS pag Fh # (1 i: Sb aint ah aren 
' i ei == \ . cay 5 . n t 
4 = 7 Zz" : 7 - | ee 7” waw 5 4 : b £ 
e a { ‘ v > _ _ a : ‘ —— ‘ —- rere anemee arate a I rn 
i . \ ne, ue 


‘tas 40 oe Ay cl i 


: =p ’ : af \ 7 . Ags 4 * 
- ; yi v 7 a 
z J ' a 


MsUG Mash Beet ass 1st hw, sae si ‘sy ibe 


— ai t ’ 
iw) Qld thai Lyin, Fh Abed rertaahed noipet 
ee) : my is : ; Preemrel — 7 i 
i ; - i} Ww i fy . 

i Ww Sieg Se a” a Gag 

: Wee [i <=tt yw W 

4 % oR a 

y ao an ® . ¥ 


Equation 21 and 22 will only be valid within a certain 


range of water saturation values. It should be clear 


from Figure 2 that She must be less than 1-Sor° 


Sw 
S 
WN Se 
0 
0) X, Xe > X 1.0 
FIGURE 2: SCHEMATIC REPRESENTATION OF WATER 


SATURATION DISTRIBUTION BEFORE BREAKTHROUGH 


The distance at which Sy reaches its maximum, 1-Sorer 


can be found from equation 20: 
we 


The foregoing treatment depends on the assumption that 
no flow occurred perpendicular to the main direction of 
fluid movement. van Meurs and van der Poel observed that 
the fingers did not deviate too much from the X-direction. 
However, Faia ie noted that especially at the top of the 


fingers there should have been components in directions 
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perpendicular to the X-axis. This phenomenon was analysed 
in more detail in the following manner. 

The formation of water protrusions at the finger 
tips requires a certain amount of water; however, once 
Pid led, these protrusions do not constitute a passageway 
for the water. This means that S|, Must exceed Som! and 
the theoretical water saturation distribution should 
therefore be cut-off at a point, Xue where the ae ren 
is larger than Som! see Figure 2. This critical water 
Saturation is denoted So: Obviously, the vertical cut- 
off as indicated in Figure 2 is an idealization. The 
real change is Seopably a gradual transition. 

To determine this cut-off point the authors reasoned 
as follows. 

The amount of water entering the unit cross section 


in one second at Xae (dis =s ' is used to £ill the form- 
Woac 


ation to the critical saturation. If the velocity at 


which S. propagates is denoted Va then: 


eV. S (24) 


(q.,) e 


Sane 


The velocity at which the critical saturation moves is 


defined by: 


) (25) 
S =S 
W (@ 


or using the definitions of X and Ws: 
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V (ee (26) 
Equations 19 and 26 combined yield: 
df 
gq | w) 
Von a Pe (27) 
o 6 cS S_=S 
Wal uC 


Substitution of Equation 27 in Equation 24 gives for the 


condition at the cut-off boundary: 


= 


qi df 

Cais: ae Ridsu (wy ae) 

q WwW S =S 
Wea 


Combination of Equations 13, 18 and 28 yields the water 
saturation at the cut-off point: 


BS 


a See 


c wa  /M-1 2» 


An expression for ae is obtained by combining Equations 
21 and 29: 
MBW 
Xo Some as ee (30) 
VE(M= 1) Sua: 
[B + /B(M 1)S. a! 
Breakthrough will occur when the critical saturation, Sar 
arrives at the outflow end, X, = 1. ‘The totalcumulative 
production at that moment, Wy is equal to the break- 


through recovery,.N.,.. Equation 30 solved: for Nop can be 


pb 
written: 


N = 28 (31) 
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where 


A= Bot VBIMET)S (32) 


Equation 31 shows that when M > ~, N +5 


pb wm?" ecre 


breakthrough occurs at the moment the minimum water satur- 
ation is established throughout the formation. 

The cumulative oil recovery after breakthrough, hae 
is easily found by integrating the water saturation 


present in the formation, see Figure 3: 


1.0 
[soe 
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FIGURE 3: SCHEMATIC REPRESENTATION OF WATER 
SATURATION DISTRIBUTION AFTER BREAKTHROUGH 
sO 
= os ax 8}6))) 
Ni, Xy (1 Sl + | S. ( 
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by substitution of Equations 22 and 23 and integration 


Equation 33 becomes: 


No =, S. + gy (2/W{MB/- W, — B) (34) 
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Equation 34 is valid up to the moment that Xy = lor 


according to Equation 23, for Ww. < MB. The recovery is 


if Ww. 2 MB. 
i€ a 


given by Np er i So 
The theory presented above was derived for a porous 
medium containing no connate water, and it involved numer- 
ous simplifying assumptions. However, in spite of these 
limitations, the theory could be used to predict the obser- 


ved behaviour of the CWO oil. The calculations and the 


comparison can be found in Chapter 7. 
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4, EXPERIMENTAL EQUIPMENT 


Viscosity correlations were obtained from the 
Research Council of Alberta for the refinery fractions, 
and the manufacturer provided this information for the 
Dow Corning fluid. The viscosities were checked with a 
Bendix Lab Viscometer. Satisfactory agreements were 
obtained. 

Interfacial tensions were measured with a Du Nouy 
tensiometer (A..Kruss, Nr. 1385). The measured values 
were corrected using the data of Zuidema and Waters~-. 
Several tests were performed on each system, and the 
reported interfacial tensions are arithmetic averages. 

Imbibition tests were conducted with apparatuses 
Similar to that described by Bobek et ala os Small 
Lucite cylinders contained the unconsolidated sand. 

The displacement equipment is shown schematically 
in Figure 4. Pictures of the experimental layout are 
presented in Figures 5 - 7. Water and oil in the bombs 
were displaced with mercury driven by one or two pos- 
itive displacement Ruska pumps. An Hg - detector con- 
nected to a relay assured that the power supply was cut 
off if mercury reached the tip of an electrode installed 


in a vertical flow line just above the bombs. This 


design removed the danger of getting mercury into the 
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BACKPRESSURE REGULATOR 


FIGURE 7 


AND SETTLING TANK 
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sandpack. 
Pressure drops were measured on Heise gauges (range 
0-300 psig). The vacuum pump was a Duo-Seal, model 1405. 
Two different coreholders were employed: 
Runt 1-46: Stainless steel tube, Length = 117.3 cm 
| ID = 5.08 cm, .OD = 6.03 cm 
Run 47-49: Lucite tube, Length = 134.8 cm 


ID = 5.08 cm,: OD =.6.35 cm 


A new type of endpiece was designed for the stainless 
Steel coreholder, see Figure 8. The Teflon packings are 
expanded when the conical metal rings are forced to move 
by tightening. the nut. The design was tested to 2000 psig 
without leaks. The Lucite coreholder has been described 
in detail by Kloepfer!?, 

Movement of the sand was prevented by a 200-mesh 
Stainless steel screen clamped down on top of the spreader 
plate, see Figure 8. The screen also kept the grooves in 
the plate free from sand so that the spreading ability 
was retained. 

The propane was driven by carbon dioxide. A back 
pressure regulator (Grove) maintained the pressure neces- 
sary to keep propane liquefied. Separation of the two 
phases was done in an old propane tank, see Figure 7, and 
the gas led to a fume hood. 


The fluid bombs were taken out of the holders and 
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filled externally. The 90° bend in the tube connections, 


see Figure 5, simplified this procedure. 
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Be EXPERIMENTAL MATERIALS 
5.1 Fluids 


Distilled water was the displacing fluid in all 
runs. It was dyed with potassium permanganate in the 
experiments with the Lucite coreholder. 

Displacement tests wate conducted with three dif- 
ferent oils. Two refinery fractions, CWO and MCT 5, were 
obtained from Imperial Oil through the Energy Conserva- 
tion Board's Laboratory. Private communication?? indicated 


that the approximate compositions of these oils in volume 


% were: 
Aromatics Naphthenes Paraffins 
CWO 25 70 5 
MET.5 15 15—80 LO= 15 


The third oil was a Dow Corning 200 minid with a 
viscosity of 100 cS. This is a silicone oil with a density 
very close to that of water at room temperatures. The vis- 
cosity is a weak function of temperature compared with 
normal hydrocarbon fluids. 

Viscosity,as a function of temperature for all_four 
fluids, is given in Table 1. The values of the interfacial 


tension between water and each of the oils are also given. 
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TABLE 1 


FLUID PROPERTIES 


Viscosity (cp) 


- Temperature 
Dow Corning 
(°F) Distilled CWO 2007 BLaid MCT 5 
Water | (100 CS) 
Es005 156 LILA 6 18) - 
0.992 WS a! 110.9 34.3 
0.980 ez 1L10%.0 33.9 
0.966 14.6 hOB sel: 33,5 
0.954 14.4 106.3 S250 
0.941 138 104.4 Sule 7) 
0.929 a 10225 30.8 
Oil Interfacial Tension Between The Oil And 
Distilled Water At 74°F 
Type (dynes/cm) 
CWO Deel 
Dow Corning 34.5 


MCT 5 34.4 
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The fluids were separated and used again. This 
should assure that equilibrium was reached in case of a 
reaction between the two phases. Viscosities and inter- 
facial tensions were measured on the recycled fluids, but 
the variations from the values obtained with fresh fluids 


were always less than the experimental errors involved. 


52°. Sand 


The coreholders were packed with a clean dry Ottawa 
sand screened to contain grains between 80 and 120 U.S. 
mesh (Fisher Scientific S-151). The narrow range of 
particle size and the subrounded nature of the grains fac- 
ilitated reproducibility of porosity and permeability. 
Uniform grain size also reduced the unavoidable inhomogen- 


eities. 


5.3 Cleaning Agents 


Propane driven by carbon dioxide was used to clean 
the sand. It was preceeded by cyclohexane to clean up 
the silicone oil, since Dow Corning reported excellent 


solubility of the 100 CS oil in this hydrocarbon liquid. 
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6. EXPERIMENTAL PROCEDURES 


Gk Packing of Coreholders 


The coreholders were packed with dry sand. Electric- 
al hammers vibrated the tubes as they were slowly filled. 
A rubber hammer was used every 15-20 centimeters for 
additional power. The vibrators were then left on for 
10-12 hours. A couple of blows with the hammer afterwards 
could produce a settling of several centimeters. The pro- 
cedure was repeated until no further settling was observed. 

The stainless steel coreholder was opened after run 


# 46. No settling or channeling of the sandpack could be 


detected when it was cleaned out. 


6.2 Initial Conditions, Porosity and Absolute Permeability 


The packed coreholder was mounted in the displace- 
ment apparatus. The Ruska pump was used to fill the inlet 
system with water up to the last valve before the sandface. 
All valves except those leading to the vacuum pump were 
closed, and this pump was left on for about 20 minutes. 
(Later experiments showed that the absolute pressure of the 
system was reduced to around 5-7 mm Hg). The valves lead- 
ing to the vacuum pump were then closed, the pump shut off 


and the inlet valve opened. The Ruska pump was used to 
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push water from the bomb into the pore space. The outlet 


valve was opened as soon as the pressure indicators on the 
pump started to climb. Water was injected until a stab- 
ilized flow pattern, input equal to output, was establish- 
ed. This usually required 2-3 pore volumes. The scale on 
the pump gave the volume injected. The amount of fluid 
ejected was measured, and the difference between the two 
was taken as the volume of the pores. The bulk volume was 
calculated and the porosity determined. 

Pressure drops across the sandpack were measured at 
several stabilized rates, and the absolute permeability to 
water calculated from Darcy's Law. The reported values are 


the arithmetic average for each case. 


6.3 Initial Water Saturation and the Effective 


Permeability. to O11 


The core, 100% saturated with water, was then flooded 
with oil at a fairly low rate (used 100 cc/hr). After the 
injection of 1000 cc, the rate was increased to 1120 cc/hr. 
and kept constant until water ceased leaving the core. 

The pressure drop was recorded and used to calculate the 
effective permeability to oil at the initial water satura- 
tion, denoted Koi- The initial water saturation was cal- 
culated from a volumetric balance. 

The core was flooded with oil after each experiment 


to restore the initial conditions. Earlier work by 
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Maguss~4 and some preliminary runs showed that this had to 


be done exactly the same way each time for any given water- 
flood series to obtain reproducibility. Even then problems 
were encountered as is evident from Figure 9. These are 
further discussed in Chapter 8. The procedure used was 

the same as the one given above to arrive at the initial 
saturation. An additional 200-1200 cc of oil was usually 
required to reach the stage where almost no water was 


produced at 1120 cc/hr. 


6.4 Waterflooding 


Rates were selected utilizing the gear facilities 
On the Ruska pumps. Production was gathered in graduated 
cylinders and centrifuge tubes. Pressure drops across the 
coreholder and the produced volumes of oil and water were 
recorded at different intervals. The breakthrough of water 
was taken as the time the first drop of water was produced 
in the effluent for the runs on the steel coreholder. At 
high rates however, some water is produced along with the 
oil almost immediately. Breakthrough was then taken as the 
point where the water-oil ratio increased by a factor of 
approximately 10. 

Runs were terminated when water-oil ratios greater 
than six were reached Because of the excessive time in- 


volved, some of the slow experiments were only brought to 
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the breakthrough point. The same was done with several 
tests performed to check the reproducibility of results. 
For each oil a residual oil saturation was obtained, and 
the effective permeability to water at that condition was 
calculated. This was the last test performed before the 


sand was cleaned in preparation for a new test series. 


6.5 Cleaning Procedures 


The sand was thoroughly cleaned before a change to 
a new 0il took place. Liquid propane was driven through 
the pack by carbon dioxide. An appropriate back pressure 
assured that the propane did not vaporize. 

Cyclohexane was employed to obtain proper cleaning 
after the runs with silicone oil. 

The complete clean-up procedures were: 
1. After CWO - runs: 3000 cc of propane at 150 psig 

3° 13000 ec. ot propane at 400 psig 
2- After Dow Corning 200 Fluid runs: 
3000 cc of ‘cyclohexane 


3 x 3000 cc of propane at 400 psig 


In both cases dried air was passed through the sand for 
three days before a new test was started. 

The reproducibility of the initial conditions, see 
Table 2, is within experimental errors. This indicates 


that the cleaning procedures were sufficient. 
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6.6 Imbibition Tests 
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Imbibition tests were carried out with small Lucite 


coreholders. 


the desired fluid. 


They were evacuated before saturation with 


The saturated cores were placed in 


the imbibition cells, and the apparatuses filled with the 


imbibing liquid. 


The amount of liquid exchange after 72 


hours was recorded. 


Tests were performed on the following 


Core System 


systems: 
1. Clean 
2. Clean 
3. Clean 
4, Clean 
oa Clean 
Fluid 


6. Clean 


sand 
sand 
sand 
sand 


sand 


sand 


and 
and 
and 
and 


and 


and 


water 
water 
water 
CWO 


Dow Corning 


MCT +5 


Imbibing Fluid 
CWO oil 
Dow Corning Fluid 
METS 
Water 


Water 


Water 
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7. RESULTS AND DISCUSSION 


7.1 Sandpack Properties 


The design of the stainless steel coreholder is 
such that the length of a sandpack has to be determined 
after packing as the total length of tube and endpieces 
less the fixed length of the endpieces. 

Length and diameter for both coreholders are given 
in Table 2 together with the porosity and the permeability 
values determined initially and after each cleaning pro- 
cedure. 

Runs 1 - 46 were all performed on the same sandpack. 
The sand was thoroughly cleaned and dried before a change 
of the displaced fluid and a redetermination of the initial 
conditions took place. The variation in porosity values 
from these tests are within the errors involved in determ- 
ining the volumes of injected and ejected water. The 
permeability variations are also within the experimental 
error. 

The porosity and permeability value for the exper- 
iments with the Lucite coreholder, runs 47 - 49, differed 
from those obtained on the steel coreholder, even though 
the same type of sand was used in both cases. However, 
the sand used to fill the Lucite coreholder was from a 


different shipment than that used for the steel tube. A 
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TABLE 2 


CORE PROPERTIES 


Diameter -= 2" {= -5.08.-cm 
all runs: 
The sand was Fisher Scientific 
S-151 80 - 120 mesh 
run a =. 46... Length = 107.38 cm (steel) 
run 47 - 49: Length = 134.78 cm (Lucite) 
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Porosity Absolute permeability Sand 
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possibility exists that the variation in sand parameters 
was caused by a distributional variation within the spec- 
ified grain size range. 

Another, and more likely cause of the higher porosity 
value was the settling problem encountered when the Lucite 
coreholder was packed. The tube was discoloured from pre- 
vious experiments and had to be washed with various hydro- 
carbon solvents. The chemicals attacked the inner surface 
and made it rough, and it was evident that this created a 
settling problem. A looser pack would explain the higher 
value obtained for the porosity with this coreholder. 

The chemicals also removed the hard glossy surface 
and made the tube vulnerable to further attack as mani- 
fested by small amounts of dissolved plastic in the 
effluent water. Plastic particles might have caused a 
‘blocking of some of the available flowpaths, and this 
could explain the lower value obtained for the absolute 
permeability. 

The porosity and permeability values given in 


Table 2 are in the same range as those determined by 


Kloepfer!?, who performed tests with the same sand. 


Pee Displacement Tests 


Pertinent data from the displacement tests are given 


in Table 3. Detailed results from each run are included in 
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Appendix B. The calculated values of the dimensionless 


group I are also given in Table 3. 
We2el — Initial» Conditions 


The initial water saturation varied between consec- 
utive runs for the three oil types tested. The variations 
encountered are plotted as a function of consecutive run 
numbers for each oil in Figure 9. 

The CWO and the MoT 5 oil can be seen to exhibit 
Similar trends, probably a reflection of their similar 
chemical composition (see Chapter 5). 

The initial water saturation was almost constant 
during the Dow Corning runs according to Figure 9. However, 
it must be mentioned. that it was necessary to sequentially 
decrease both the amount of oil injected and the injection 
rate to prevent a decrease in the value of Si’ The water 
repellency of the silicone oil was believed to be the 
cause of this phenomenon. The imbibition tests performed, 
see Appendix A, showed that sand saturated with Dow Corning 
fluid imbibed less water than the two other sand-oil 
systems. Water had difficulty in entering the sand once 
it had been contacted by the silicone oil. 


The problems with the initial water saturation are 


further discussed in Chapter 8. 


» 
¢ 
ie 
, 
eT | 
; 
yi d 
“8 
, “ 
“~ 
yr BS 


Pa Wh 
20 


; ty sire 


an ay ah pine , 


A ebre eel a snes anal 


i 
t 
i 
gt 
2 *a9 n 
a7 
, 
w 
“ja 
> « 
‘7 7 
= 
Bab , 
= 


ae i 
yd Re A) ek? LEO: 
’ i i iP ie 


Sead sods Udita “eit nant eh: So wee i 


Dd Th, Tig orate Ls Odd | Xe, 


AM, 


AW 


- ar pists te | pei hetee 


ye sj 


1 


tad 


Ww baveruVes bros ene? bh Bh a “Abad we? 
By > «/, ae) ad } ‘ | : 
ie dee ews aed * opel iy fe Khe cae: ersten tik 
‘ m~ i; 1B by oe Y a f ae } j is er = ab he saab a laa ar 
} ; « | j sah 
P20 OF, Cee Sey ie? Eg ‘any weer: bart: ab ay 
bap 75.) Lett tg) gale sinptart oar 
he \ 
ab ‘sia ni Sener 7p hl aot 
- ' a 4 i : eo, = 
i ~ ; n 
4 f i. ' a ee 
' i * ’ re : vy ; 


47 


M0 HOV3S YO4 SYSSWAN NNY SAILNOASNOD 
O02 G| Ol S 


Vv v Vv 
Vv 


(vl = ator ) ONINYOD MOG Mf 


(sgl = advod) OM) a @ 


(zpe=mnyol) ¢ 19; 


NOILVYNLVS Y3LVM TWILIN] SHL JO NOILVINVA 3HL ©6 3uNDIS 


(Ad %) NOILVYNLVS Y3SLVM IVILINI 


7 


ar 
Pubokd 


f° 


nee 


~ = 
— 


a, 


} 
‘3 


Ser 
“kh 


sda 


aa 
ayia 


le 


ie aay 


\ DR ee ; fie = 


r ~ 
‘= 
of 
ih j ed y 
Cee | ae: 
od ) Zo an i. fie en 
| 7 et eae 
a i" if 
a - 
iJ = d 
A’ ee. ri 
5 UY hy ty i A =| 
; ' = 
" a 
i { 
. if ay y = 
4 
Wi fs 7 i 
S a” my 40 
Cs a 
OT ie eet 3 t 
‘ ae 
SS 1 ip | _ 


48 


Note that the effective permeability to oil at Sa 


(Koj in Table 3) stayed almost constant for each oil, even 
though the initial water saturation increased significant- 
ly during the CWO and the MCT 5 runs. This can be taken 

as support for the reasoning presented in Chapter 8 on the 


nature of the initial water saturation and its effect on 


the performance of the tests. 


ospeiera Ok Recovery at Breakthrough, The Effect of 


Emulsification and Viscous Fingering 


Oil recovery as a function of cumulative amount of 
water injected at low, intermediate and high rates for 
each of the oilsemployed is presented in Figures 10 - 12. 
‘The breakthrough points are marked with arrows. The cor- 
rection applied to the breakthrough point for the Dow 
Corning curve in Figure 12 is commented upon later. 

The lines representing recovery before breakthrough 
have varying slopes in these Figures due to the differences 
in the initial water saturations. 

The recovery at breakthrough decreased from the low 
rate to the intermediate rate for all three oils studied, 
and it was significantly higher for the oil with the lowest . 
viscosity. The breakthrough recoveries for the MCT 5 and 


the Dow Corning are about the same for the tests shown. 
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The relative importance of the production after 
breakthrough increases with increasing rates and viscosity 
ratios according to the CWO and MCT 5 curves in Figure 10 
and 11. The Dow Corning fluid shows the same variation 
with rate, but the effect of the higher viscosity ratio 
is not evident when compared with the two refinery oils. 

Figure 12 shows a different relationship between 
the two refinery fractions. Recovery at breakthrough is 
still higher for the CWO oil, but more of the MCT 5 was 
recovered at any value of cumulative water injected. This 
improvement in the recovery of the MCT 5 oil was observed 
in all runs where the displacement rate was higher than 
200 cc/hr (i.e... runs 335, 38,. 43, 44 and 46). These runs 
correspond to values of the scaling group I larger than 
0.2 in Figure 14. It is evident from the sudden increase 
in the breakthrough recovery value observed in Figure 14, 
that a change in the mechanisms controlling the displace- 
ment took place for a value of I between 0.2 and 0.4. 

The change is believed to have been caused by the 
formation of water-in-oil emulsions in the pore structure. 
A single phase "slurry" was produced after breakthrough 
in all these experiments. The stability of the emulsions 
varied, and the time required for the separation of the 


two phases was anywhere from 10 minutes to 10 - 15 hours. 
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The increased breakthrough recoveries can be explain- 
ed as an interaction between two factors which affect the 
recovery: 

1. The viscosity ratio between the displacing front 
and the displaced oil decreases since the vis- 
cosity of an emulsion depends mainly on the vis- 
cosity of the external phase”, 1.e. the oil -in 
this case. 

2. The emulsification taking part at the displac- 
ing front leads to a partial blocking of paths 
with low resistance to flow, and this can be 
expected to improve the microscopic sweep 
efficiency. 

The data shown in Figure 14 indicates that no furth- 
er improvement of breakthrough recoveries could be obtained 
once an I value of 0.8 had been reached. The breakthrough 
Fecoveryiatsthic point, 3/7. 7%) Om t.O 4. boy 4s seen aconbe 
approximately the same as that obtained with a very low rate. 

Figures 13 to 15 show the recovery at breakthrough 
versus the log of the scaling factor I for the three oils 
employed. 

The data from the CWO and the MCT 5 experiments 
have been split: into two groups representing runs perform- 
ed in the early stage, where the initial water saturation 


still changed rapidly between each experiment, and the 
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late runs performed with each oil, where the change in 

Sy eas less dramatic. The split can more or less be 
anticipated from Figure 9. In the case of the CWO oil the 
first 8 runs are seen to have an initial water saturation 
less than 18.5%, where as the last 12 runs have S,, values 
between 20.4% and 25.6%. Only the latter points were 
considered when the curve in Figure 13 was drawn. The 
first 8 runs with the MCT 5 had Sad values less than 24%, 
whereas Sd varied between 27.7% and 33.5% for the last 

9 runs which were utilized for the drawing of the curve 

in Figure 14. 

The data obtained with the Dow Corning fluid, 
Figure 15, could be correlated without discarding runs. 
This could be expected in view of the almost constant 
value for Soi during these experiments (see Figure 9). 
However, it was necessary to correct the breakthrough 
value obtained in run 29. Breakthrough determined as a 
10 fold jump in WOR values gave an abnormally low value 
for this run. The reason for this was probably the 
previously mentioned problem of retaining the Soi level. 
The water in the system became easier to move with every 
run, and this resulted in a premature breakthrough de- 
termination for run 29. The recovery history for this 
experiment was shown in Figure 12. The breakthrough 


value given in Table 3 is indicated together with a 
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a corrected breakthrough value obtained by considering an 
extension of the straight line relationship up to this 
point. Both points are plotted in Figure 15. 

The curves from Figures 13 - 15 have been trans- 
ferred to Figure 16 to facilitate comparisons. Figure 
17 shows breakthrough data from the work of de Haan”, 
Engelberts and Klinkenberg~ and Collins?! plotted versus 
the log of the same scaling number. 

All the curves shown in Figures 16 and 17 exhibit 
the same features if the anomally of the MCT 5 curve, the 
emulsification effect at I values larger than 0.2, is 
ignored. Similar curves were also obtained by Jones- 
Parra et Eye and Newcombe et a for comparable 
systems. 

Low values of the scaling group I yielded high 
breakthrough recoveries. This is explained by the in- 
stability theory as the region where the peak-to-peak dis- 
tance is larger than or in the same order of magnitude 
as the tube diameter. The displacing front will therefore 
always move in a piston-like fashion and high breakthrough 
recoveries can be expected. 

Breakthrough recovery decreased with increasing I 
values in an intermediate range. This is attributed to 
the formation of viscous fingers. Oil was bypassed and 


left as "islands" in the pore space behind the flood front. 
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The breakthrough recovery stabilized for I values 
larger than a fixed number for each system. This means 
that the amount Oe Oil bypassed by the front remained 
approximately constant as the fingers became more numer- 
ous and the distance between the peaks decreased. 

Values for the constant C in Chuoke et al.'s* 


theory can be calculated for the systems shown in Figure 


16. Recall Equation 4 (Chapter 2): 


fee (ee hey) De 
= Ow 
c= ey os a (4) 


and the dimensionless group I determined by: 


aT: = soe 
yvK 


a combination of these two equations yields: 


31 i /it BPD: 
LYK 


(35) 


The value Of I in Equation 35 must: be taken, at ithe 
point where the decrease in the recovery at breakthrough 
becomes eee: A point was assigned to each curve as 
indicated by the arrows in Figure 16. The C values cal- 
culated are given in Table 4, together with the value 


used by de Haan? to determine the points shown in Figure 
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TABLE 4 
SYSTEM ViSCH -RARTIO © 
CwO a ewe) Ziel 
MCP'5 34.2 24.3 
Dow Corning La A 39.4 
de Haan 4 & 25 30 


The values of C obtained for the three systems 
employed are seen to be very close to that used by de 
Haan? to determine the point where the formation of vis- 
cous fingers started to influence his breakthrough recov- 
ery data. He took the value from Chuoke et ales” 
Original paper where they found that a C value of 30 gave 
the best fit between the theory and their visual obser- 
vations. de Haan's paper and this study indicate that 
Equation 35 can be used to estimate the value of I at 


which the viscous fingering phenomenon becomes significant 


by assuming C = 30 for a wide range of viscosity ratios. 


7.2.3 Recovery After Breakthrough 


Baugure! 18 isda plot of ‘the recovery obtained atua 


water-oil ratio of approximately six (6) versus the log 


Ot ss 


The data are seen to be scattered. The lines for 
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the CWO and the MCT 5 oil were drawn by discarding the 
early runs as done before. Only two data points were 
then left for the MCT 5 case, so this line was dashed. 
The Dow Corning data points were not exhibiting any 
trend. 

The recovery at WOR = 6 is seen to increase with 
increasing values of I for the CWO experiments. The 
dashed MCT 5 line indicates the same, but suffers from 
the lack of further data points. 

Figure 19 shows the recovery at a cumulative 
throughput of water of 0.5, 0.3 and 0.4 pore. volumes for 
the CWO, the MCT 5 and the Dow Corning oil respectively 
olettea asa function of the log of the scaling number I. 

It appears as if the recovery obtained at an equal 
amount of injected water stabilized for each oil some 


time after breakthrough. 


7.2.4 Applying the Theory of van Meurs and van der Poel” 


The following equations were derived in Chapter 3: 


Recovery at breakthrough was given by: 
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and the recovery after breakthrough was determined by: 


2 1 MB - te 
Se gay (200 ae) B) (34) 


A value for = fae was obtained by injecting water at 
a high rate until practically no more oil was produced. 
This was done after the last run with each oil, and the 


values obtained can be found in Appendix A. 


a) The Cwo System 


From Appendix A: Sor = 0.204 


Runs where viscous fingering occurred were selected, 


and the value of S taken as the S_. value in Table 3. 
wm wi 

ihe Viscosity ratio, M, was takenvas 15.5 70m all runs. 

The calculation of the breakthrough recovery, Ne is 

performed in Table 5 and values obtained experimentally 


have been entered for comparison. 


TABLE 5 


CALCULATION OF RECOVERY AT BREAKTHROUGH 
FOR THE CWO OIL 


Run hit B A oe Exp. (Pe Vv.) 
10 0.211 0.5385 Pages 0.408 Mere ei? 
Le On228 9 0.568 L938 O2427 0%.336 
PZ 0.217 0.579 1.929 OVS 0.366 


ute) 03233 0.563 Wao ae 0.432 Yes) 
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The calculated values are seen to depend heavily on 
the values introduced for Som’ van Meurs and van der Poel 
found that see = 0.15 gave the best fit for the data avail- 
able to them. That would yield Nes = 0.335 for the CWO 
experiments if Sor is assumed to be the same as above. A 
value of S| between 0.15 and 0.211 leads to theoretical 
breakthroughs between 33.5% and 40.8% of the pore volume. 
This encompases all the breakthrough recoveries obtained 
in the viscous fingering region with this oil. 

‘The theoretical recovery after breakthrough was 


calculated from Equation 34 for run 10, since this showed 


the best fit at the breakthrough point: 


Run 20 (320 cc/hr) 


Ws 0.451 0.648 0.874 
EXp. Ny 0.403 0.458 0.482 
Catc. Ny 0.426 0.468 O5.05 
A a5 0.023 0.010 0.023 


The theory predicts the observed behaviour surpris- 
ingly well considering the many assumptions involved in 
arriving at Equations 31 and 34. 

However, the theory could not be used to predict 
the behaviour observed with the Dow Corning and the MCT 5 


oil as is evident from the following. 
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b) The Dow Corning System 


From Appendix A: SA = 0.47 
From Table 3: S ZAEOe LO 
wm 
From Table 1: Mo Suda el 
this yields: B = 0.43 
A = 2.619 
and Nob = 04143 P.v. 


The experimental values varied from 16.9 to 29.6% P.V. 


according toetabile 3: 


c) MmMcT 5 - System 


From Appendix A: See! 0.147 
Eromsvabley 3: i S420 528 
From Table 1: M = 34.2 
this yields: B = OF 573 
A = 2.915 
and N = 0.434 P.V. 
pb 


The breakthroughs obtained experimentally were be- 


tween 14.1 and 25.9% P.V. 


7.2.5 Results From the Transparent Lucite Core, 


Runs 47 - 49 


The Lucite coreholder and the dyed injection water 


were employed to be able to visually observe the behaviour 


‘ x I irene’? Batseyv: Poet, ee | 


hy 
, eno : a 
tn eee 
ee Sh oe 
: a0e 6. = 
J HER. i 
j wee LA 
ef af 8, SS Steck ahead 9 na ite adwvonstrdaesd ait aS 


ee cor epee ai eee bas Lik 


4 
cs SS At) i Jt peg D abe, 
- aay Sa > 
1 
‘i 7) mn ¥5 = . 
a TPO DISSE bio! ev 496% 1S 
. He 


Bel Gs Vez? BD Seiko? yi trike By ida 2 3, BB 


of the CWO system. Run 47 was performed with clear water 


to. stabilize the fluid saturations. 
a) Run 48," Rate” =» 6.25 -cc/hr. 


No fingering was observed at this rate. Some clear 
water was produced just ahead of the arrival of the dyed 
water at the outlet end. Approximately 10 cc of oil was 
produced between the arrival of the dyed water at the end 
of the sandpack and the first appearance of coloured water 
in the effluent. 

This indicated that the error introduced by the out- 


let effect was around 1.3% of the pore volume for the 


determination of breakthrough with the CWO oil at low rates. 


b). Run -# 49, Rate. = 320 ¢ce/hrs 


| 

Numerous fingers were formed at the early stage of 
the flood. However, they broke up and formed a transition 
zone with a tongue along the bottom part of the tube. The 
behaviour was similar to that described by Perkins and 
Johnston*t, 

Clear water was again produced, but this time over 
the whole period before breakthrough. 

Only about 2-3 cc of additional oil was produced be- 


tween the arrival of the dyed water at the outlet end and 


its appearance in the effluent. This indicated that the 
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error introduced by the outlet effect was in the order 

of 0.4% P.V. for breakthrough determinations at this rate. 
The dyed water was successfully washed out after the 

completion of run 48. Approximately 2 pore volumes of 

clear water was used before the sand was resaturated with 


oil in preparation for the next run. 


7.2.6 Applying Buckley - Leverett Theory 


The possibility of applying the Buckley-Leverett 
theory to obtain the recovery at breakthrough was suggest- 


ed. The average mobility at this point can be obtained 


from: 
K L 
Ce). = (36) 
Rie ee 
where 
S) = average mobility of the fluids at 
breakthrough 
AP, 1 = pressure drop across the core at 
breakthrough 


It was also mentioned that Equation 36 could be 
used as a criterion to determine the recovery at break- 
through. Theoretically experiments are only comparable 
aide (K/u) is the same for each run. The calculations 


are performed below for runs 10 - 13: 
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Run g Sa (K/U) 
(cc/hr) (psia) 

10 320 ae as9 

11 480 8.6 Lissa: 

LZ 640 ls: ees 

3 960 3:8 eS: 


Note that the following conversion is necessary to 


obtain consistent units above: 


ae A 
(=), = se NA Sh 19021642, —= 
60° x 20.268 AP, om. AP, om. 


These results indicate that the scattering in the 
data presented might have been caused by the inability 
to obtain the same conditions for each run as evident from 
the random changes in the average mobilities calculated 


above. 


8. ERROR SOURCES 


8.1 Experimental Errors 


The mercury in the pump cylinder was subjected to 
the static head of the mercury reservoir when filled and 
the static head of the fluid bomb after each injection 
cycle. The volume of injected fluid taken from the scale 
on the pump will therefore be in error due to the com- 
pressibility of the mercury. The error was found to be 
in the order of 1-4 cc for the porosity determinations, 
depending on the levels of mercury. This would show up 
as a maximum error in the porosity values obtained for 
runs 1 - 46 of approximately 0.2%. 

The respective volumes of oil and water in the pores 
at any moment were thereafter calculated from volumetric 
balances. The total volume ejected was measured and taken 
as the true volume injected. The possibility of a cumula- 
tion in the error discussed above was thereby prevented. 

Gravity segregation affected these tests since the 
oils had densities different from that of the water. The 
effect was reduced by rotating the coreholder 180° each 
time the injection fluid was changed. However, the core 
could not be rotated during an experiment. The tests 
performed on the transparent system indicated that gravity 


segregation might be an important factor even though it 
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has been neglected in most of the work done on horizontal 
linear displacements. 

The permeabilities calculated are subjected to errors 
involved in determining the pressure across the sandpack. 
The pressure drops during the determination of absolute 
permeabilities were in the range of 2 - 3 psia, and small 
errors in zero-setting or reading off the gauge could 
easily explain the variation from 16.7 é6 17.1 darcys 
obtained in the tests on the steel coreholder (see Table 
2, Chapter 7). 

The determination of breakthrough for runs 1 - 46 
was influenced by the outlet end effect. The magnitude 
of the error introduced was in the order of 0.4 - 1.3% P.V. 


according to the experiments with the transparent system. 


8.2 The Nature of the Initial Water 


Run 48 and 49, conducted with dyed injection water, 
showed that some of the initial water in the system was 
mobile. Water was always produced along with the oil at 
high rates (see the tables in Appendix B), but at low 
rates no water was produced until just before breakthrough. 
This difference in behaviour was probably caused by the 
capillary end effect, which prevented water from leaving 


the core until the exerted pressure gradient exceeded the 


capillary pressure. 


The mobile nature of the connate water has been 
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SES It has been shown 


studied by several researchers 
that connate water, in many cases, is banked up ahead of 
the displacing fluid, and it is this bank which actually 
displaces the 0i17°'29, 

Breakthrough determined as the point when water 
reaches the outlet end differs from the actual breakthrough 
of injected water when the initial water is banked up 
ahead of the front. It was impossible ro estimate to what 
degree this affected the breakthrough recoveries obtained. 

The recovery at breakthrough was seen to increase 
with increases in the initial water saturations during the 
preliminary runs with both the CWO and the MCT 5 oil when 


_the curves in Figures 13 and 14 were established. Two 


explanations of this behaviour seems to be possible: 


1. A contact period was necessary to bring out 
effects of the water-wet nature of the sand. It 
retained more and more water closely bounded 
to the grain surface as time passed, and less 
water was left to move and cause premature 
breakthroughs. 

2. W/o emulsions were formed in the sandpack. The 
water globules in the oil phase would be dif- 
ficult to move, and the imerease inethes initial 
water saturation will therefore not lead toa 


greater mobility of the initial water. 
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Emulsions could have been formed with the Cwo 
oil in view of the similarity in composition 
between the two oils. However, these must have 
been unstable since no emulsification was evi- 


dent in the effluent during runs 1 - 20. 


An equilibrium stage would be expected if reactions 
between the sand grains and the fluids were important. It 
is evident from Figure 9 that the Si never reached an 
equilibrium level in the CWO and the MCT 5 runs. However, 
it had a tendency to stabilize for a few runs at various 
levels. 

It might be that the variations observed could be 
explained by a combination of 1 and 2 mentioned above. 

There is no doubt about the fact that the greatest 
uncertainties in the data obtained in this study were 
introduced by the variation of the initial water satura- 
tions between consecutive runs and the mobile nature of 


this water. 
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9. CONCLUSIONS 


The following conclusions are made based upon the 


experimental data and the discussion presented: 


ls 


The initial water in an artificial unconsolid- 
ated core is mobile. The determination of 
breakthrough as the moment water reaches the 
outlet end might therefore differ from the 
value obtained if the breakthrough of the inject- 
ed water could be detected separately. 
Formation of viscous fingers caused decreasing 
breakthrough recoveries. This took place for 
values oe the scaling factor I between 0.005 
and 0.02 with the three oils studied. 

High breakthrough recoveries were obtained at 
low values of "I". 

The breakthrough recovery stabilized at high 
values of the scaling factor. The stabilization 
took place for "I" values larger than 0.1 with 
the MCT 5 and the Dow Corning systems. "I" had 
to exceed 0.8 in the case of the CWO oil before 
breakthrough recoveries stabilized. 

The value of the scaling group I at the point 
where viscous fingers will begin to form can be 
estimated for all three systems studied if a 


value of 30 is assigned to the constant C in 
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Chuoke et Bales theory. The same value was 
obtained by them in their experimental work, and 
de Haan” used it to successfully predict the 
onset of fingering in his experiments. 
Breakthrough recoveries were drastically im- 
proved when a water-in-oil emulsion was formed. 
The simplified theory of van Meurs and van der 
Poel might deserve more attention. The problem 
seems to be the determination of reasonable val- 


ues for S and. So. 
wm or 


17 


28W Site ates aa 


p RIGS < . ; i+; mace J Se i ay Prd ote: ¥ 
RT 1 Pep S iia & 4 Fu tillenagiey 


tka 


et ree 


J ee \ 
) oI Sie Tete eee 


4 ee ae ; 

“3 ( aes hy 
tA oF Pea Sy Sy) ey eae f 
nea a 


SLO, Sh AGN Oe tenia se) we 


ee See ’ a) 
= Vf = As - i" 
’ beg i 
; } Piet vat 
oh t by ‘ a9 i i > , ‘is i 
; t m PL R y oe i Ps , pee ae " j " a * , 
ay 7 ! a rt. = ¥ 7 teoteti ‘3 at at 
oe i + ; ; i : Po, oe oe 5 
Ter - > u a 
* 2 sy iy ‘a or) ‘4 aw 16 : y 
: ae - | a 
Y ( . Mi 7 oth ‘ a A l ; a 
5 1 ! Py Pees By a“ 
= a OE ‘7s Bs A ae ye 


10. RECOMMENDATIONS 


The following modifications of the displacement 

apparatus might prove useful: 

1. Design a high pressure swivel connection so 
that the coreholder can be rotated during a 
displacement experiment. A slow rotational 
movement would reduce the observed effect of 
gravity. 

2. The outlet end of the stainless steel core- 
holder should be redesigned, if possible, so 
that breakthrough can be detected as the 
moment the injected fluid reaches the end of 
the sandpack. 

The following suggestions are made for further 

experimental work in the same area: 

3. Conduct an investigation of the qature of the 
initial water. The effect of the mobile char- 
acter on the determination of breakthroughs 
should be studied, together with the effect of 
letting cores "age" before they are employed 
in a displacement experiment. 

4. Consider the necessity and feasibility of dying 
either the injected water or the water used to 


establish the initial saturation. 
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Tests should be performed with systems having 
viscosity ratios in the order of 500 and higher 
in view of the current interest in the recovery 


of extremely viscous crude oils. 


The diameter of the coreholder should be varied 
to study the effect of the canal width on the 


occurrence of viscous fingers. 
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dimensionless constant 
diameter of tube 


fractional flow of water 


LVu 

W ; 
——— , dimensionless group 
yVK 


absolute permeability, om? 


absolute permeability, darcy 
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shock mobility ratio 
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pressure in oil phase 
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PROPERTIES OF FLUIDS AND SAND 


The Gasoline and Oil Laboratory, Research Council 
of Alberta measured the following viscosities and densities 


for the CWO and the MCT 5 oils: 


Viscosity In Centistokes 


cwo McT 5 
at 70°F 12:69 39°53 
at 100°F 9.47 19.04 
at 130°F 593 10.74 

Densities 

cwo MCT 5 
at 60°F eat sy cit: 0.8618 
at 70°F 0.8581 0.8581 
at 100°F 0.8474 8006474 
© API Sy eg 


The viscosity data was plotted on a standard viscos- 
ity-temperature chart (A.S.T.M. D 341-43) to obtain the 
kinematic viscosity, n, at room temperatures. The den- 
sities, op, were plotted vs temperature on normal graph 
paper, and the dynamic viscosity at each temperature cal- 


culated from: 
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The values are entered in Table 1, Chapter 5. 

Viscosities for the Dow Corning 200 Fluid were cal- 
culated from the specification charts provided by the 
manufacturer. 

The viscosity of distilled water was obtained from 
Perry's "Handbook of Chemical Engineering". 

The results of the imbibition tests after 72 hours 


were as follows: 


System Imbibed Volume (cc) 
1. Sand & CWO imbibing water 0.3 
2. Sand & D-C ee water Ord 
3. Sand & MCT 5 imbibing water OnS 
4, Sand & water imbibing CWO 0.0 
5. Sand & water imbibing D-C 0:60 
6. Sand & water imbibing MCT 5 0.0 


These results indicate that the three systems were water- 
wet. However, the sand saturated with Dow Corning 200 
Fluid, imbibed only one-third of the water amount observed 
with the two refinery fractions. The water repellancy of 


the silicone oil prevented water from entering to the same 


degree. 
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The residual oil saturation, See and the permeabil- 


ity to water at this point, Kig 1 were determined after 


or 


the completion of the last experiment with each oil. The 


following values were obtained: 


Oil-type Sor (3) Kis (darcy) 
OF 
After run 20 CWO 20.4 4.16 
After run 29 Dow Corning 47.0 2.28 


After run 46 MCT 35 14.7 Soul b 
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